نشریه 1 سال بیست و پنجم» شماره یک, ۱۳۹۲ 
کاربرد الگوی ۷۷۸۲ در شبیه‌سازی عددی امواج غلتشی * 


)0 ۱ بر 0 
علی مهدوی ناصر طالب بیدختی" 


چکیده در اين مقاله الکوی شار متوسط وزن‌دار (۷۸۳) برای حل معادلات غیر حط یآب کم‌عمق و بررسی تکامل و انتشار امواج غلتشی 
د رآبراهه‌های با شیب تند بهکار رفته است. در حقیقت این الکو دفت شارهای محاسبه‌شده توسط روش کلاسیک سر را افزایش می‌دهد. 
همچنین الگوی "۷۷۸ علاوه بر توانایی شبیه‌سازی ناپیوستگی‌های مرتبط با امواج غلتشی, یک رهیافت قوی برای از بین بردن ناپایداری‌های 
عددی ناشی از عمق ناچی زآب, که معمولا در شبیه‌سازی امواج غلتشی مشاهده می‌شوند, اراله می‌نماید. نیم‌رخ‌های محاسبه‌شده عم قآب و 
سرعت جریان تطایفی رضایت‌بحش را با حل تحلیلی موجود نشان می‌دهند. نتایج عددی مربوط به تکامل زمانی دامنه‌ی موج برای جریان- 
هایی با اعداد فرود مختلف با آنچه که از یک مدل دیکر به‌دست م یآید, مفایسه شده‌اند. هم‌چنین مفایسه‌هایی میان نتایج ایین دو مدل برای 
تاریخچه‌ی زمانی عم قآب صورت پذیرفته است. تطابق مشاهده شده, از کارایی و دقت الگوی پیشنهادی در مدل‌سازی امواج غلتشی حکایت 


می‌کند و این در حالیست که الکوی حاضر در پیاده‌سازی رایانه‌ای ساده‌تر است و به زمان اجرای ی کوناه‌تر نیاز دارد. 


واژه‌های کلیدی معادلات غیرخطی آب کم عمق. موج غلتشی. الگوی تسخیر شوک. شار متوسط وزن‌دار. 
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۰ ۲۷۸۱ 5۱۵607 6۱0۵۳۵۵۲ 6۵6۵6 000 ۱۵6۵۵00 موم ( عآمرورآه 


۵۵6 ۲7۵۱۵64 ,مدمه مصتساویه مه ۷27۵ ۲۵1 رعممتامناوه مه ۲«مالمطو تعممتامماد . ق16۵۲ 
۱۹ 


۷ تاریخ دریافت مقاله ٩۰/۷/۲۵‏ و تاریخ پذیرش آن ٩۱/۱۰/۲۶‏ می‌باشد. 
(۱) نویسنده‌ی مسوژول. دانشجوی دکتری. بخش مهندسی را ساختمان و محبط زیست. دانشکده مهندسی. دانشگاه شیراز. 
(۲) استاد بخش مهندسی رام ساختمان و محبط زیست. دانشکده مهندسی. دانشگاه شیراز. 
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با در نظر گرفتن جریانی یکنواخت بر روی بستر دارای 
اصطکاک در یک مجرای باز شیب‌دار. هنگامی که عدد 
فرود این جریان دست‌نخورده‌ی یکنواخت از مقدار 
مشخصی بیش‌تر شود ناپایداری‌هایی در سطح آزاد 
جریان رخ می‌دهد. تحت چنین شرایطی چنانچه یک 
اغتشاش کوچک به جریان دائم یکنواخت اعمال گردد. 
در نهایت مجموعه‌ای از امواج شکننده با آشترک 
(00۲0) در الگویی پله‌مانند در جریان به‌وجود می‌آید. 
این الگوی ناپیوسته‌ی متناوب به امواج غلتشی ( ]۲01 
0۳۹ موسوم است. امواج غلتشی با سرعت ثابت و 
تغییرات متوالی رژیم جریان منتشر می‌شوند. جریان 
فوق بحرانی با گذر از یک پرش هیدرولیکی در یک 
چارچوب متحرک با سرعت یکنواخت به جریان زیر 
بحرانی تبدیل می‌شود (شکل ۱). به بیان دیگر از 
دیدگاه ناظری متحرک که با سرعت موج غلتشی به 
سمت پایین‌دست جریان حرکت می‌نماید. نیمرخ موج 
با تغییرات متوالی جریان فوق بحرانی به زیر بحرانی 
همراه است. این در حالیست که از دیدگاه یک ناظر 
ثابت. جریان در سراسر آبراهه‌ی فوق بحرانی می‌باشد. 
امواج غلتشی به‌صورت معمول در مجاری ساخته- 
ی دست انسان نظیر آبراهه‌ها و سرریزها به‌وجود 
می‌آیند. این امواج در فلوم‌های آزمایشگاهی نیز تولید 
شده‌اند [1]. رخداد چنین امواجی علامت این است که 
تغییرات در جریان و عمق آب می‌تواند قابل توجه 
باشد. هم‌چنین وجود امواج غلتشی ممکن است سبب 
سرریزی از کناره‌های آبراهه گردد. این امر مشکلاتی 
عملی برای مهندسان هیدرولیک به‌همراه دارد [2]. 
هرچند امواج غلتشی اغلب در مجاری مصنوعی رخ 
می‌دهند این امواج در جریان‌های طبیعی نیز دیده 
شده‌اند که از آن جمله می‌توان به آبراهه‌های یخی 
جریان‌های ثقلی در آزمایشگاه» اقیانوس و دریاچه‌ها 


اشاره نمود [3]. علاوه بر این اغتشاشات شبیه به امواج 
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غلتشی در موارد گوناگون دیگری نظیر جریان چند 
فازی [4], جریان گل [5] و جریان در مجاری کشسان 
مانند جریان هوا و خون در رگ [6]» رخ می‌دهند. 

بخش عمده‌ای از تحقیقات مرتبط با امواج غلتشی 
به تعیین شرایط لازم برای شکل گیری امواج غلتشی 
اختتسامن استه سسکا وهافته تسایر انز فده 
عبارت است از بررسی پایداری یک جریان یکنواخت 
بز رزوی شیب بت که اعتطاشنی بنز سسطع رذآ 
اعمال گردیده است. در مطالعات پیشین, باحل 
معادلات ناویر-استوکس برای جریان لایه‌ای و اعمال 
اغتشاشات سینوسی بر سطح آزاد جریان. عدد فرود 
آتنتانه برای رخداد موج غلتشی در حدود 20.5 ,۲ 
تفن له اشت :81 رز هه ای تعر نان ها 
آشفته, با در نظر گرفتن یک آبراهه‌ی عریض, توزیع 
سرعت یکنواخت و ضریب اصطکاک ثابت. عدد فرود 
آستاندق رخداد موج غلتشی برابر با 22 مآبه دست 
ای [ 9 

باید توجه داشت که رخداد امواج غلتشی به وجود 
اصطکاک در مقابل جریان بستگی دارد و در صورت 
عدم وجود اصطکاک موح غلتشی شکل نمی گیرد 
هم زان شراب رم بای تک رایخ 
امواج را بر مبنای مقاومت جریان به‌دست آورد [11]. 
با این وجود امواج غلتشی بر روی بسترهای بسیار 
نامنظم به‌وجود نمی‌آیند[12]. به بیان دیگر مقاوست 
پیش از حد جریان مانع از شکل گیری امواج غلتشی 
می‌گردد به طوری‌که برای تشکیل امواج غلتشی بایستی 
مقاومت جریان از یک مقدار مشخص کم‌تر باشد 
[13]. شایان ذکر است پدیده‌ی امواج غلتشی می‌تواند 
مانسیل فرسانش خاکهرا افانش فد دلیل انتن امتر 
تغییرات زمانی و مکانی جریان در خلال رخداد چنین 
امواجی گزارش شده است [14,15]. 


سال بیست و پنجم. شماره یک ۱۳۹۳ 


ما ۲۱۱۷۱۵۲۵۵۱۱۵ 
من : 


٩0‏ اسمه 


1۵ 


سل 


شکل ۱ انتشار امواج غلتشی در یک آبراهه که با تغییرات متوالی جریان از حالت فوق بحرانی(5100) به حالت زیر بحرانی(570) 


در یک چارچوب متحرک از طریق ایجاد پرش هیدرولیکی همراه است [20] 


به احتمال قوی 1۳01128 اولین محققی بود که 
امواج غلتشی بزرگ دامنه (06انامصقءع »1 را 
به‌صورت تحلیلی توصیف نمود[11,16]. وی با در 
نظر گرفتن یک موج متناوب با شکل و سرعت ثابت. 
به تولید نظری نیم‌رخی شبیه به نیم‌رخ امواج غلتشی 
پرداخت. بر مبنای ایده‌ی بنيادین این محقق, 127655167 
حلی بسته برای امواج غلتشی دائمی در یک آبراهه‌ی 
مرش با امیظگای کات ازاته کان | تا سس ۶ 
آن 97001 به توسعه‌ی یک مدل نظری بر مبنای 
معادلات غیر خطی آب کم عمق ( 202110627 
فصمتامناوه ۳/2/6۲ 9۳21100) پرداخت [17]. این 
معادلات که بر فرض توزیع فشار هیدرو استاتیک و 
سرعت متوسط گیری شده در عمق استوارند برای 
شبیه‌سازی جریان‌هایی که در آن بعد عمودی جریان از 
بعد افقی آن بسیار کوچک‌تر است. به کار می‌روند. 
تاکنون محققان بسیاری بر پایه‌ی این معادلات به 
بررسی جنبه‌های مختلف شکل‌گیری و انتشار امواج 
غلتشی پرداخته‌اند. 
#0 و همکاران (۲۰۰۵) به کمک معادلات غیرخطی 
آب کم‌عمق به بررسی جریان شیاری (1107 ) 
پرداختند [18]. این محققان بر مبنای کارهای گذشته 
[13,17] و با فرض برابری عمق متوسط جریان و 
عمق جریان یکنواعت» ضمن ارائه‌ی روابطی برای 
نیمرخ و سرعت موج نشان دادند که پدیدار شدن 
امواج غلتشی, افزایش تنش برشی جریان را به‌همراه 
دارد و از همین روی سبب تسریع فرآیند فرسایش 
خاک و تولید شیارهایی در سطح خاک می‌گردد. 
(2007) 00070062 در چارچوب معادلات غیر خحطی 
آب کم‌عمق. شکل گیری امواج غلتشی در سیلاب های 
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جریان یافته بر بستر شیب‌دار را مطالعه نمود [19] 
این محقق با به‌کارگیری روش تفاضل محدود برای 
حل معادلات حاکم» تکامل خطی و غیرخطی اغتشاش 
اولیه را بررسی نمود. (2006) 26 20 0۵6 با ارائه‌ی 
یک مدل حجم محدود که در آن از الگوی موسوم به 
(016ظ) (رعاممیک صرح وومعن) رتهع2008ظ) برای 
حل معادلات آب کم‌عمق استفاده شده بود. به بررسی 
شکل‌گیری و تکامل امواج غلتشی پرداختند [20]. با 
وجود توانایی الگوی 13016 در شبیه‌سازی جریان‌های 
مختلف. این الگوی عددی در پیاده‌سازی رایانه‌ای 
پیچیده است و در مقایسه با الگوی مورد استفاده در 
این مطالعه به زمان اجرای طولانی‌تری نیاز دارد. 

در این مقاله پس از تشسریح معادلات حاکم 
جزئیات روش عددی مورد استفاده برای تخمین شار 
عددی و چگونگی پیاده‌سازی عبارت چشمه در مدل 
عددی به تفصیل ارائه می‌گردد. در ادامه پس از بیان 
معیار ناپایداری جریان و شرایط اولیه و مرزی برای 
شبیه‌سازی موج غلتشی, با طرح یک آزمون کلاسیک 
شکست سد و مقایسه‌ی نتایج با حل تحلیلی موجود. 
توانایی مدل در شبیه‌سازی یک جریان ناپیوسته مورد 
ارزیابی قرار می‌گیرد. هم‌چنین مقایسه‌ای از نتایج مدل 
عددی و حل تحلیلی برای یک موج غلتشی ارائه 
خواهد گردید. آن‌چه کم‌تر در متون علمی پیشین مورد 
توجه قرار گرفته‌است چگونگی شکلگیری نیم‌رخ 
سطح آزاد در حالات آستانه‌ی رخحداد و عدم رخداد 
موج غلتشی می‌باشد. به‌همین منظور این پدیده‌ها در 
انتهای این نوشتار بررسی خواهند شد. 
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معادلات حاکم 


معادلات آب کم‌عمق دارای کاربردهای وسیعی در 
مهندسی هیدورلیک و مهندسی سواحل می‌باشد. این 
معادلات که در واقع قوانین بقای جرم و اندازه حرکت 
را بیان می‌کنند. تاکنون در شبیه‌سازی پدیده‌هایی نظیر 
امواج جذر و مد در مصب‌ها و آب‌های ساحلی, انتشار 
آتراخ رکه تفر بااروی آسرام بای پر 
هیدرولیکی ایستا و جریان رودخانه و مجاری باز, با 
موفقیت به‌کار گرفته شده‌اند. معادلات آب کم‌عمق در 
حالت پایستار (10110 0601561۷211۷70) و یک‌بعدی به 
شکل زیر نوشته می‌شوند: 


۱( (۷) بر 2۷ 


0) 0 

رابطه‌ی فوق در واقع دستگاهی از معادلات 
غیرخطی هذلولوی بوده که در آن ۷ بردار متغیرهای 
پایستان (7)۷ بردار شار (۷۵۵۱۵۲ ال و ٩0۷(‏ 


عبارت چشمه (61170] 5011706) نامیده می‌شوند: 


1 9 
<(۳0۷ ,| | ۷ 
تب ۳0( ۳ 
)۲( 


«#9 ۳ ِ ۷ 


در رابطه‌ی فوق ع شتاب گرانش.» مختصات 
افقی» ] زمان, 1-96 عمق آب» ٩,‏ شیب بستشس 
(اباحنا سرعت افقی متوسط گیری شده در عمق و 
0۶ قاریتضط خاک مت مس باشد. [20] ممادلتش 
(۱) تحت عنوان صورت دیفرانسیلی معادلات بقای 
جرم و بقای اندازه‌ی حرکت شناخته می‌شود و استفاده 
از آن تنها در مواردی که جواب‌ها هموار می‌باشند 
صحیح است. هنگامی که یک ناپیوستگی در جریان 
وجود دارد» بایستی از حالت انتگرالی معادلات استفاده 
شود. معادله‌ی (۱) پس از انتگرال گیری در حجم کنترل 

به‌صورت زیر قابل بیان است [21]: 
طسیرط ‏ - ۳ 


1+1/2 1-2 +۹ )۳۲( 
۸ ۸ : 
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که در آن ‏ ۵اندازه‌ی سلول محاسباتی و ۸۵۶ گام 
زمانی می‌باشد. پایین نوشت 1 معرف مرکز سلول 
محاسباتی 1 (:<) و بالانوشت‌های ۲ و 11 
به‌ترتیب نشان‌دهنده‌ی گام‌های زمانی معلوم و مجهول 
ورربزگ و مر نیز که به شارهای بین سلولی 


موسومند, به ترتیب به شار عددی در نقاط ورر,:< > 


(مرز بین سلول های 1 و 1+1) و ورر< (مرز 
بین سلول‌های 1-1 و ز) اشاره می‌نمایند. در قسمت 
تعلربه چگونگی تخمین این شارهای عددی پرداخته 


می‌شود. 
الگوی شار متوسط وزن‌دار برای معادلات آب 


کم‌عمن 

تاکنون برای حل عددی دستگاه معادلات پایستار 
روش‌های گوناگونی طراحی شده و به‌کار رفته است. 
این روش‌ها عموماً در قالب روش‌های اجزای محدود 
المان محدود و حجم محدود طبقه‌بندی می‌شوند. از 
این میان روش حجم محدود از آن‌جایی‌که در حبل 
جریان‌های ناپیوسته به عبارت‌های اضافی نظیر لزجت 
مصنوعی نیاز ندارد و حاصیت پایستاری معادلات را 
به خوبی حفظ می‌نماید از مقبولیت بیشتری برخوردار 
است. 

روش‌های گودانف (306000058 00000۷) که با 
حل مسأله‌ی ریمان در مرز هر سلول همراه می‌باشند و 
از جهات انتشار موج برای تخمین شار بین سلولی بهره 
می‌برند. یکی از مهم‌ترین زیرمجموعه‌های روش حجم 
محدود به‌شمار می‌روند. 37001 و همکاران (۱۹۹۹) 
ضمن استفاده از روش گودانف برای حل جریان در 
لوله‌های تغییرشکل پذیر نشان دادند نتایج حاصل در 
مقایسه با الگوی تفاضل محدود مک کورمک 
(0۷20)00102016) از دقت بیش‌تری برخوردار است 
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در اين مطالعه برای تخمین شار عددی از الگوی 
شار متوسط وزن‌دار (۷۷۵۳) ( ۸۷۵۲۵۵6 ۱۷۷۵1۵0160 
[۳) که روشی از نوع گودانف است و یک حل‌کننده 
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تقریبی مسأله‌ی ریمان به‌شمار می‌رود. استفاده گردیده 
رای ای ی تن وی در موه بان 
به‌صورت ترکیبی وزن‌دار از شارهای مرکزی سلول‌های 
چپ و راست مرز و شار در ناحیه‌ی میانی موسوم به 
ناحیه‌ی ستاره تخمین زده می‌شود «شکل۲). شار 
عددی در ناحیه‌ی ستاره ("/۳)3» با استفاده از حل 
کننده‌ی تقریبی ریمان از نوع رال ر عضا رجها:۴2 
۲ ۷۵ 200) با توجه به شارهای چپ و راست 


مرز و سرعت امواج به‌صورت زير محاسبه می‌شود. 
- "0۷ 
٩۲) ۷۱-۹۳۲) ۹۵۵‏ 
- 5 


() 


که در آن ,۷-۶۷ و ,,۷< ۷ هم‌چنین ٩,‏ و 
ی ی 
مرز 1+1/2 اشاره می‌کنند. برای تخمین این سرعت‌ها 
از روابط زیر استفاده می‌شود [23]: 


("طول+ ای ول دهع ٩‏ 


( دول ای ول اطع ٩‏ 
)0( 
در عبارات فوق "1 و "لا به ترتیب مقادیر عمق و 
سرعت جریان در ناحیه ستاره می‌باشند که به صورت 
زير قابل محاسبه‌اند: 


2 
1 1 | 1 ب 
0 ی ۳ 
1 
ول - نو + (ودا + 3 ۷ 
)0( 


برای ساختاری شبیه به آن‌چه در شکل (۲) نشان 
داده شده است. شار ۳ در مرز 11/2 به صورت 
11/2 11/2 


بل 
9 رت ۳ 
121 


ال بیست و پنجم شماره یک ۱۳۹۲ 


2۷ 


2/۸ ,۰ که در آن ٩,‏ سرعت موج 6 


1 
)۸( اک و . ات 8 ما ات 2 1 


در روابط (۸ و۷ تعداد معادلات پایستار 
(برای دستگاه معادلات یک بعدی آب کم‌عمق: 
2 ۲) و ورب بیانگر شار عددی در بازه‌ی 6 به 
طول ۷,۸ است (شکل ۲). با اعمال شرط موسوم به 
تغییرات کلی از بین رونده (1۷۲) ( عمنامزته۷ 10121 
مصنطمتمتصز» که متضمن جوابی عاری از نوسان غیر 
فیزیکی در نزدیکی نواحی ناپیوسته جربان است. 
رابطه‌ی (6 به شکل زیر بازنویسی می‌گرده [21]. 
ی 
اه 


۳۹۵ 
مش 0( )ول < - 
9 


(8) 

که در آن ( )520 به تابع علامت اشاره می‌نماید. 
با تعریف شارهای (2۳)۷ ور ( 2۳6۷ و3 و 
(۳)۷ 5 ۳ پرش در شار عددی به شکل 
وررب۳- ۲ < و۵۳ محاسبه می‌شود. در رابطه 
(» وراه تابع محدود کننده شار جهت اعمال شرط 
می‌باشد که با اعمال مقدار مناسب استهلاک 
عددی به الگوی ۰۷۷۸۳ یکنوایی (۵050101010) 
حل در نواحی دارای گرادیان زیاد را تضمین می‌نماید. 
تاکنون انتخاب‌های گوناگونی برای تابع محدود کننده‌ی 
شار پيشنهاد شده است [21]. در مطالعه‌ی حاضر از 
محدود کننده‌ی موسوم به ۳۳اه استفاده شده 
اتکنت* 

(به|-1)-1< 0 


۱۰ 
9" ات 2) 0 2 ,)9 ۱ 
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شکل ۲ ساختار الگوی ۷۷۸۳ در مرز سلول محاسباتی برای حل معادلات یک‌بعدی آب کم‌عمق 


که در آن: 


۸۳ 
< 6 101 اس 
۸ 
) ۱ 0 1+12 لت 
(۳ط۸ 


13/2 


00 


بل تن ۱ 
22 < 0) پرضص 


3 )( 

در رابطه‌ی فوق, ,+ 7 

در متغیر * در خلال موج 6 می‌باشد [21]. شایان 
ذکر است که روش ۸۵۳ با اعمال اندکی تغییر در 
سرعت‌های موح. به‌راحتی برای مسائلی که بخشی از 
بحشک می‌باشد قابل استفاده است. به‌عنوان نمونه 
می‌توان از این روش برای شبیه‌سازی مرز متحرک خحط 
ساحلی و مطالعه‌ی پدیده‌های انتشار و بالاروی امواج 


پیاده‌سازی عبارت جشمه در الگوی عددی 
آن‌چه در بخش پیشین به تفصیل بیان گردید. تنها برای 
0 < (50۷ ( قابل استفاده است. در حضور عبارت 
چشمه می‌توان الگوی ۷۷۸۳ را بدون تغییر به‌کار برد 
به‌شرطی که عبارت چشمه طی یک یا چند گام 
جداگانه انتگرال گیری, وارد محاسبات گکردد 
[24 ,21]. در مطالعه‌ی حاضر برای دست‌یابی به‌دقت 


مرتبه‌ی دوم در زمان» دستگاه معادلات زير که متشکل 
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از سه معادله‌ی مقدار اولیه است در الگویی متوالی حل 


ک قم 
۱7 
۱ ۷ - -- :01۳5 
239 
۷ :105 
۳۲ 02۷ 
0-<--+_- :۳[(۳6۶ 
۲ ۵ب اش 2 
. ا. :105 
۱7 
2-۷ -- :0۵۳5 
۶ ۱ (۷) 3 
. :105 


که در آن ۸۲۸۰/2 می‌باشد. فرآیند حل 
سه‌مرحله‌ای دستگاه فوق بدین نحو قابل تشریح است: 
در مرحله‌ی اول با استفاده از مقادیر جریان در گام 
زمانی قبل به‌عنوان شرط اولیه. یک معادله‌ی دیفرانسیل 
معمولی که در برگیرنده‌ی عبارت چشمه است. حل 
می‌گردد و شرط اولیه الگوی ۷۷۵ به‌دست می‌آید. 
خروجی الگوی ۱۷۸۳ در مرحله‌ی دوم به‌عنوان شرط 
اولیه در مرحله سوم برای حل معادله‌ای نظیر مرحله‌ی 
اول اعمال می‌گردد تا مقدار متغیرهای جریان در گام 
جدید زمانی حاصل شود. به‌طور حلاصه این فرآیند 
حل به‌شکل رهام (اش)پر( شام - اور قابل 
بان اب کا ور او و 1 یرس عم کر عتساوت 
چشمه و عملگر معادله‌ی همگن نامیده می‌شوند. در 
مطالعه‌ی کنونی از الگوی انتگرال گیری ذوزنقه‌ای برای 
حل معادلات چشمه (مراحل اول و سوم) استفاده 
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کردنده استد این عملگر ضمنی دارای دقتی از مرتبه‌ی 
دو است و برای معادل‌ی مرحله‌ی اول با رابطه‌ی زیر 


داده می‌شود (پیوست الف): 


۳ (۹0۴ ۵ - ه ۶ 
از 2 2۱۱ 


که در آن ] ماتریس واحد است و ۷۵-۷۳ - ۸۷ 
معرف پرش در متغیرهای پایستار می‌باشد. عبارت 


ماتریس ژاکوبین بردار چشمه می‌باشد. 


معیار پایداری محاسبات 
از آنجایی که الگوی ۷۸۳ الگویی صریح می‌باشد» 
پایداری محاسبات با اعمال معیار کورانت- فردریش- 
لویی(,۳1)) بر مقدار گام زمانی تأمین می‌شود. 


ره 1 ,06 >0 صنص ,) < ۸۱ 


۸ 
رطل برس 


که و ان بل) عدد بآتل) می‌باشد. رابطه‌ی فوق 
بیش از طول یک گام مکانی را طی کند. 


شرایط هیدرولیکی جریان ناپایدار 

حل جریان یکنواخت به‌صورت << و مدا < ا 
برای معادلات غیر خطی آب کم‌عمق از توازن میان 
نیروهای گرانش و اصطکاک بر روی یک شیب ثابت 
حاصل می‌شود به نحوی که: 
)6 شردام < رگر(ت 

امواج غلتشی معمولاً بر روی آبراهه‌های دارای 
شیب تند باعدد فرود جریان بزرگ‌تر از ۲ 
( ,2/20 < ودا) اتفاق می‌افتند [1]. این موضوع با 
توجه به رابطه‌ی (۱۵) به‌شکل زیر قابل بیان است. 
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1۹ 


٩, < 4, ٩ 
رابطه فوق که نشان دهنده‌ی شرایط ناپایدار جریان‎ 

است در حقیقت بیانگر عدول از معیار پایداری جریان 
یکنواخت ( مطعل2 ک ولا) می‌باشد. این معیاز براساس 


شرایط اولیه و مرزی برای شبیه سازی امواج 

همان‌طور که قبلاً اشاره شد امواج غلتشی با اعمال یک 
اغتشاش به جریان ناپایدار به‌وجود می‌آیند. با در نظر 
گرفتن شکل سینوسی برای این اغتشاش, شرایط 
اولیه‌ی عمق جریان عبارت است از: 
۱۷( | (,)صنوع +1] وط ع (0,)ظ 

که در آن , به حاصل‌ضرب عدد مسوج 
(۵۷60000۲) در عدد 27 اشاره می‌نماید و ع 
عامل بزرگنمایی (100107 21001111621100) اغتشاش 
او یه ی پاش که ترافس شین اعتاق 
به عمق جریان اولیه‌ی یکنواخت. در محاسبات صورت 
پذیرفته در این تحقیق عامل بزرگ‌نمایی برابر با ۸۰/۵ 
ین اوله‌ی,جریان یکتواعت در نظر گرفیه فده ات 
به‌عبارت دیگر 0.005 2 می‌باشد. هم‌چنین سرعت 
اولیه‌ی جریان به شکل زیر اعمال می‌شود [20]: 


)۱۸( ( ,0+ کیک)ط ۵ + وا ع (۷0),۵ 

که در آن دو متغیر ,9 اختلاف فاز میان عمق و 
سرعت اغتشاش و ضریب , هر دو توابعی از سرعت 
اولیه‌ی جریان یکنواخت و عدد مسوج و فرکانس 
زاویه‌ای اغتشاش می‌باشند [20]. در مدل عددی. 
شرایط مرزی در دو انتهای چپ و راست دامنه‌ی 
محاسباتی به‌صورت دوره‌ای (06110016) در نظر 
گرفته شده است. بر این اساس در هر لحظه مقادیر 
متغیرهای جریان در دو مرز انتهایی با یکدیگر برابر 
می‌باشد. به‌عبارت دیگر: 
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۱2 ۳ 


که در آن ,۶ و لا به‌ترتیب به مختصات مرزهای 
چپ و راست ناحیه‌ی حل اشاره می‌نمایند. در تحقیقی 
(2002) :اتعحاهما مه طعتاانطع2 ب‌کمک روش 
۸ به مطالعه‌ی امواج غلتشی پرداختند [25]. مدل 
این محققان به داده‌های سرعت و دامنه‌ی موج غلتشی 
به‌عنوان شرط مرزی بالادست نیازمند می‌باشد و این 
درحالی است که مدل توسعه يافته در مطالعه‌ی حاضر 
با توجه به شرط مرزی اعمالی رابطه‌ی )۱1٩(‏ به 
داده‌های یادشده نیازی ندارد و ایين شرط مرزی با 


خاصیت دوره‌ای امواج غلتشی هم‌خوانی دارد. 


شکست سد. بررسی امواج ناشی از شکست سد. یک 
مسأله‌ی بسیار جالب نظری و در عین حال کاربردی در 
مهندسی هیدرولیک می‌باشد. به‌عنوان اولین آزمون 
محاسباتی الگوی عددی ارائه‌شده برای شبیه‌سازی 
امواج ناشی از شکست سد به‌کار می‌رود. در این مثال» 
آب ساکن بالادست به عمق 5 21.0 ,9 به‌وسیله‌ی 
یک دیواره‌ی قائم (سد) که در وسط کانالی به طول 
0ب( گرار عارته از ان اکن بایین دست به 
عمق ۲ 20.05 ,9 جدا می‌شود. در این حالت نسبت 
عمق آب پایین دست به بالادست برابر با 
3ن رظارظ هی تاش و مرزهای سخبدو واست 
دامنه‌ی محاسباتی با شرط مرزی غبر بازتابی ( 208 
0 در مدل عددی گنجان ده می‌شوند. در 
لحظه‌ی 1-0 دیواره به‌صورت ناگهانی برداشته 
می‌شود و یک موج آشترک به سمت پایین‌دست و یک 
موج منفی به سمت بالادست انتشار می‌یابند. مقایسه‌ی 
نتایج حاصل از شبیه‌سازی عددی برای نیم‌رخ سطح 
آزاد جریان در زمان 8 21.0 با حل دقیق [21] در 
شکل (۳- الف) نشان می‌دهد که الگوی عددی ارائه 
شده به خوبی قادر به شبیه‌سازی ناپیوستگی عمق در 
جبهه‌ی آشسترک می‌باشد. در شکل (۲- ب) روند 
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افزايش دقت مدل با زیاد شدن تعداد سلو‌های 
محاسباتی به‌روشنی قابل مشاهده است. هم‌چنین نتایج 
شبیه‌سازی سرعت جربان در زمان 5 21.0 در 
شکل(۳- ج) با جواب تحلیلی مقایسه شده است. 
تطابق میان نتایج این آزمون و حل تحلیلی, تأکیدی بر 
ویژگی تسخیر شوک مدل ارائه‌شده می‌باشد. 


وااوج 100 س س - 
وااهه 200 ۰-۰- ۱ 1 
جااوج 400 سس 
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شکل ۳ الف: نیم‌رخ طولی موج ناشی از شکست سد. ب: 
جزئیات جبهه‌ی موج ج: سرعت جریان در زمان 1.05 -] 
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شبیه‌سازی موج غلتشی 


در این بخش الگوی عددی پیشنهادی برای مدل‌سازی 
استفاده قرار می‌گیرد. یک جریان دائنم یکنواخت با 
شرایط 2.5<و "وه 0.001ع<مو و 0۸<0.006 در 
لحظه‌ی 1-0 تحت اغتشاشی با مشخصات 2107 ,1 
و 0.005 < و قرار دارد. شیب بستر بنا بر رابطه‌ی (۱۵ 
برابر با 0.0375 2 ,5 محاسبه گردیده و دامنه‌ی 
محاسباتی به طول 22.050 با اختصاص ۲,2۶0 
و 2.01 < ۲ تعریف می‌شود. در شکل(4) تکامل این 
مختلف به تصویر کشیده شده است. پس از گذشت 
زمان موج به تکامل نهایی خود رسیده و قطار 
موج متشکل از پنج موج متوالی بدون تغییر شکل با 
سرعت ثابست در طول آبراهه به پایین دتییتت 
انتشار می‌یابد. نیم‌رخ این قطار موج و سرعت جریان 
در شکل(۵) با حل تحلیلی اراثه‌شده در [13] 
مقایسه شده‌اند. بیشینه خحطای محاسباتی (تعتی : 


8 1 


آمعتارامحه هه ظ 


و )۷) در حضیض موج 
رخ می‌دهد و برابر با 10/۷ می‌باشد. 

همان‌طور که در بخش مقدمه بیان شد. شکل گیری 
امواج غلتشی با تغییرات متوالی رژیم جریان در طول 
آپراهه همراه است. تحت همان شرایط جریان و 
اغشاتن پنشین» چنانچه ناظری, با مدرعنت نوج غلاشتی 
با جریان حرکت نماید. عدد فرود جریان از دیدگاه وی 
مطابق با آنچه که در شکل (0) به تصویر کشیده شده 
تغییر می‌نماید. تحت شرایط مذکون سرعت ابت 
حرکت موج غلتشی برابر با "عط 20.55 می‌باشد. 
هم‌چنین نتایج نشان می‌دهد که با افزایش عدد فرود 
نان وله و ومد وان که از ضرق 
ول /(0-۷)< برظ در محدوده‌ی وسیع‌تری تغیسر 
ی تمابد رک ۷ دلیل آیرن افو ایرد بت کته افیر این 
سبب افزایش ,۳ فوق بحرانی ( (08) ,در 
شکل ۷) و نیروی مخصوص (10۲66 50661116) قبل از 


پرش هیدرولیکی می‌گردد و این عامل به کاهش عمق 


سا بیست و پنجم شماره یک ۱۳۹۲ 


0۱ 


جریان توب رای فلز برتن میاه هس 
منظور بایستی نیروی مخصوص بعد از پرش نیز 
افزایش یابد تا تعادل نیروی مخصوص (با احتساب 
نیروی اصطکاک بستر) برقرار گردد. به‌روشنی افزایش 
نیروی مخصوص بعد از پرش سبب می‌شود که 
جریان زیر بحرانی عمق بیشتر و عدد ,؟ کمتری را 
تجربه نماید ( (18) ,در شکل ۷). همان‌گونه که 
نتایج نشان می‌دهد برای حالت آستانه رخداد موج 
کته بضر ان ا تن گاه تا تشر مب انس 
21 مج ۴ < (صنص) ,۲) خواهد بود. 

مقایسه‌ی تاریخچه‌ی زمانی عمق و سرعت جریان 
در وسط آبراهه (0 1.0-<) در شکل(۸) حاکی از 
تطایق رضایت‌بخش نتایج مدل حاضر با نتایج مدل 
6 220 ۵6 [20] می‌باشد. هم‌چنین اين محققان با 
فرض تغبیرات نمایی دامنه‌ی اغتشاش, یک رابطه‌ی 
لگاریتمی برای رشد و زوال دامنه‌ی موج به‌دست 
آوردند که در محدوده‌ی خطی صادق است و لگاریتم 
طبیعی دامنه‌ی موج را به‌صورت تابعی خطی از زمان 
بان با ام کفوو با کلفیی نساز و امس 
یافتن اثرات غیرخطی. رشد دامنه‌ی موج کند می‌گردد 
و در نهایت متوقف می‌شود. این موضوع در شکل(٩)‏ 
که لگاریتم طبیعی دامنه‌ی موج بر حسب زمان برای 
اعداد فرود مختلف را نشان می‌دهد. به‌حوبی نمایان 
است. شرایط جریان و اغتشاش درست مانند حالت 
قلح تفه کنر وان شا ای ترا 
زمانی داسته‌ی موج تطایقی قابل قبول با زابطه‌ی تحلیلی 
ارائه شده در [20] نشان می‌دهد. حال آن‌که با گذشت 
زمان دامنه موج به مقدار ثابتی می‌رسد. بازه‌ی زمانی 
بروز رفتار خطی موح با کاهش عدد فرود جریان 
یکنواخت. طولانی‌تر می‌شود. برای فراهم آوردن امکان 
مقایسه در ناحیه‌ی غیر حطی: نتایج مدل 3016 نیز در 
شکل )٩(‏ آورده شده است. همان‌طور که مشاهده 
می‌شود هر دو مدل در پیش‌بینی دامنه موج تکامل‌یافته 
برای اعداد فرود مختلف جواب‌های یکسانی به‌دست 


داده‌اند. 
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شکل ۱۱ مقایسه‌ی نتایج مدل حاضر و الگوی 3016 برای 
تاریخچه‌ی زمانی عمق جریان در وسط آبراهه. (الف): 
رد-۳ 


در حالت حدی ۱2۶2 که آستانه‌ی رخداد موج 
غلتشی به‌شمار می‌رود. همان‌گونه که انتظار می‌رود 
هیچ‌گونه تغییری در دامنه‌ی موج شبیه‌سازی شده 
به‌وجود نیامده است. هم‌چنین برای اعداد فرود 
کوچک‌تر از اين مقدار حدی» تغییرات دامنه‌ی موج 
نسبت به زمان دارای سیری نزولی می‌باشد. این زوال 
دامنه به همراه آستانه‌ی رخداد موج غلتشی در شکل 
(۱۰) در قالب نیم‌رخ‌های سطح آزاد جریان به‌تصویر 
کشیده شده‌ان د. به‌روشنی در حالت 21.5 ,لآ 
همان‌طور که 
دامنه‌ی اغتشاش با گذشت زمان مستهلک گشته و از 
بین رفته است (شکل ۱۰-الف»). حال آن‌که برای 


22 ,"۰ هرجند تقارن اولیه‌ی موح تاحدودی 


دست خوش تغییر گردیده است (اين موضوع در نتایج 
مدل 13016 نیز دیده می‌شود). نیم‌رخی با دامنه‌ی ثاببت 
نسبت به زمان مشاهده می‌شود (شکل ۱۰- ب). برای 
حالات نشان داده شده در شکل (۰) تغییرات زمانی 
عمق جریان در وسط آبراهه با آنچه که مدل ۳801 
به‌دست داده است. در شکل (۱۱) مقایسه گردیده 


است. برای حالت 21.5 ,3 این دو الگو زوال دامنه‌ی 
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اغتشاش را باروندی مشابه و تطابقی قابل قبول 
پیش‌بینی کرده‌اند (شکل ۱۱- الف). دو خط افقی بالا و 
پایین در شکل (۱۱- ب) بیشینه و کمینه‌ی عمق را 
مطابق با دیدگاه نظری نشان می‌دهند. بر این اساس. 
بیشینه و کمینه‌ی عمق برای آستانه‌ی رخداد موج 
غلتشی به‌ترتیب باروابط مط(ع +201 ,ررظط و 
ط( )1-6‏ ,روط قابل بیان می‌باشند. همان‌طور که در 
شکل (۱۱- ب) مشاهده می‌شود. مدل حاضر عدم 
رشد و عدم زوال دامنه (دامنه ثابست) برای آستانه‌ی 
رخداد موج غلتشی را تا حدودی به‌تر شبیه‌سازی 
نموده است. 

شکل (۱۲) لگاریتم طبیعی دامنه‌ی موج بر حسب 
پارامتر ,1 (حاصل‌ضرب عدد موج در عدد 27) برای 
اعداد فرود مختلف را نشان می‌دهد. نتایج حاکی از آن 
است که به‌ازای هر عدد فرود ثابت. افزایش ,6 با 
کاهش دامنه‌ی موج همراه است. شایان ذکر است که 
افزایش عدد موج به‌معنی کاهش طول موج می باشد. 
علاوه بر این به‌ازای هر ,6 ثابت. افزایش عدد فرود 


سبب افزایش دامنه‌ی وی غلتشی گردیده تخت : 


12 14 16 


نشریه مهندسی عمرآن فردوسی 


10 


کاربرد الکوی ۷۷۸ در شبیه‌سازی عددی امواج غلتشی.. 


(2006) جع" لصج ۵60۷-1۷]60۷7100۷2۵>نار] 
ضمن انجام مقایسه بین الگوهای مختلف حل عددی 
معادلات آب کم‌عمق بر اهمیت زمان مورد نیاز برای 
انجام محاسبات رایانه‌ای (زمان پردازنده‌ی مرکزی: 
عصسنا تا) تأکید و از آن به‌عنوان معیاری جهعت 
ارزیابی بهینه بودن مدل عددی یاد کردند [26]. 
به‌همین منظور در این مطالعه مقایسه‌ای میان زمان 
پردازنده‌ی مرکزی مدل حاضر و مدل 13616 تحت 
شرایط مشابه جریان صورت پذیرفت. برای انجام این 
مقایسه ۵۰ ثانیه از جریان با استفاده از ۱۰۰۰ سلول 
محاسباتی و عدد بآ برابر با 0.65 ,0 توسط این 
دو مدل بر روی رایانه‌ای با پردازنده‌ی مرکزی 
و وه کی با هاش کت یه 
یکسان در محیط 90 ۳01۵ اجرا شدند. 
همان‌طور که نتایج این مقایسه در شکل (۱۳) نشان 
می‌دهد. الگوی پیشنهادی در این مطالعه به زمان 
اجرایی حدوداً یک‌سوم زمان اجرای مدل 13016 برای 
حل مسائل مشابه احتیاج دارد. 


او بش 


۱0/۳ 


4 6 8 
27۳-۲ 9 


شکل ۱۲ لگاریتم طبیعی دامنه‌ی موج به‌عنوان تابعی از ,>1 
برای اعداد فرود مختلف 


سال بیست و پنجم. شماره یک. ۱۳۹۳ 


۱006۱ ۳۲۵56۲ هس 


[[ 


کت 
9 
ه‌ 


۳ 


9 
ه‌ 


۱۷۵06۱ )۵ هب | 


زره 
‌ 


(و) ۱۴0۶ لاطن 
۳ 
۵ 


شکل ۱۳ زمان مصرفی پردازنده‌ی مرکزی جهت شبیه‌سازی پنجاه انیه از جریان برای اعداد مختلف فرود اولیه 
با استفاده از مدل حاضر و الگوی ۳1 


جمع‌بندی و نتیجه گیری 

مدل غاد ارائه شنده در این مطالعه از الگوی شار 
متوسط وزن‌دار برای تخمین شار عددی و یک الگوی 
ضمنی برای پیاده‌سازی عبارت چشمه در معادلات 
یک‌بعدی آب کم‌عمق در حالت پایستار بهره می‌برد. 
مدل توسعه‌یافته بدون نیاز به هرگونه عبارت اضافی 
نظیر لزجت مصنوعی. که تعیین آن معمولاً با سعی و 
خطا توأم است. قادر است که هر گونه ناپیوستگی در 
متغیرهای جریان را به‌نحوی موثر در تعداد نسبتاً کمی 
سلول محاسباتی حل نماید. طرح یک مسأله‌ی 
کلاسیک شکست سد و مقایسه‌ی نتایج مدل با حل 
تحلیلی موجود نشان داد که علاوه بر سرعت جریان؛ 
هم جبهه قائم موج آشترک و هم نواحی ملاییم نیم‌رخ 
موج منفی به خوبی شبیه‌سازی گردیده‌اند. بررسی 
فرآیند تکامل اغتشاش و شکل‌گیری موج غلتشی در 
چارچوب یک مسأله مقدار اولیه و مقایسه نتایج حاصل 
با حل تحلیلی موجود. حاکی از توانایی الگوی 
پیشنهادی در مدل سازی عمق و سرعت جریان می- 
باشد. شبیه‌سازی جریان برای مقدار آستانه‌ی رخحداد 
ناپایداری» به ایجاد نیم‌رخی بدون رشد در دامنه 
اتجامیه نز این دی خالی استه کشیرای اعفاد کترود 
کوچک‌تر از مقدار آستانه اغتشاش به‌طور کامل از بین 
رفته و جریان پایداری خود را حفظ می‌کند. 

گر کی ترا کار یدانق مرح 


سال بیست و پنجم. شماره یک ۱۳۹۲ 


نسبت به پارامتر ,1 (حاصل‌ضرب عدد موج در عدد 
7 بررسی گردید و نشان داده شد که به‌ازای هر عدد 
افزایش می‌یابد. علاوه بر آن نتایج تعتاکی از ان اشتتن 
که با افزایش عدد فرود جریان اولیه مت عدد فرود از 
دید گاه ناظر متحرک ۳ در محدوده‌ای ریت 3و تس 
تغییر می‌نماید که اين به‌معنی عمیق‌تر شدن حضیض 
اخیر در متون علمی پیشین کم‌تر مورد اشاره قرار گرفته 
است. 

استفاده از اگوی حاضر و الگوی 801 برای 
بررسی تغییرات زمانی دامنه‌ی موج اک از تطابق 
این در حالی‌است که الگوی پیشنهادی در پیاده‌سازی 
رایانه‌ای بسیار ساده‌تر است و برای شبیه‌سازی یک 
جریان حاوی ناپیوستگی نظیر امواج غلتشی به زمان 
محاسباتی نسبتاً ناچیزی نیاز دارد. 


سپاسگزاری 
بدین‌وسیله نویسندگان مقاله از آقای ۷0-116 6 و 
پروفسور 6 سک از دانشگاه علم و فناوری هونگ 
کنگ (۳161791) به خاطر در اختیار گذاشتن مدل 


عددی ۳71 قدردانی می‌نمایند. 
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1 


ی ی 
چشمه رابطه‌ی (۱۳) تشریح می‌گردد. این عملگر 
ضمنی بسرای حسل معادله‌ی دیفرانسیل معمولی 
(5)۷ - 0۷/4۱ با شرط اولیه "2۷ ۷ به‌کار می‌رود. 
برای نیل به این هدف در اولین گام عبارت چشمه 
به کمک بسط تیلور تقریب زده می‌شود: 


تن ۴( + ۸ (فگ). (5)۷- 50۷7 


از ظرفیخ به‌کمک قاعده‌ی مشتق زنجیره‌ای می‌توان 
نوشت: 


050۷(  050۷( ۷ 

را ب ‏ 
در رابطه‌ی فوق 08)۷(/2۷. ماتریس ژاکوبین 
بردار چشمه است که با رابطه‌ی (الفب-۲) محاسبه 


می‌شود: 
0 
(الف -۲) نام 
1 


0500 _ 
0۷ 2٩ 


مشتق زمانی بردار متغیر یایستار با استفاده از روش 
تفاضل محدود قابل محاسبه است: 


فر چم 


کاربرد الکوی ۷۷۸ در شبیه‌سازی عددی امواج غلتشی.. 


۵۷ ۷-۲ 

2 
با معرفی ۷-۷۳ - ۸۷ و جایگذاری روابط 
۱۱| 


(الف-) 


عبارت چشمه به شکل زير بازنویسی می‌شود: 


(الف-۵) .۰ 00۵+ ۸۷ |22 ]+(0ه -(()/5۱ 


در همین حال طبق قاعده‌ی انتگرال گیری ذوزنقه‌ای 
(الف-۲) [(۵+(۷ 3 + ین ۷ ۷ 
با جای‌گذاری رابطه‌ی (الف-۵) در رابطه‌ی 
(الف -1) ور پس از اند کی عملیات جبری» عملگر 
چشمه برای معادله‌ی اول از دستگاه معادلات (۱۲) 
حاصل می‌گردد: 
زر 0۷ 2 


که در آن 1 ماتریس واحد می‌باشد. 


14017 .عم ,(95)4 ,7 .ایاه0 .۲ ,"عامصصمط هون صا مصتقج مبجنسلامد ۵۴ تصعصم‌ماه1۳۵۷" بیط .3 بکله20ظ .1 


(1998) ,9۵ نا ر5ع۳۲ و۸ ,۲۱۵۷۲۲ امصصف مهو ۵۶ ومتاناهت وا رگ رفعاطم]۷ . 
2004(۰) و1-33 .00 ,514 ۱۷۲66 0ذبا۲ .ل ۷۵۷۵۵ 0۴۲۵۱۱ ومتصمصنا رو تمه مره رز بل روط . 


7 1۱ طمتاقز۵؟ باه ۵۶ صوتصهطمع/۲(۷ 


.(2000) ,977-99 .00 ,26 ,1۳۱0۷ ۱۵56 م۷ . .10 ر 
۰ و(110)3 .۱۱0 .0۲۵۵1 .ل ,110۷7۲ 4ماهتاصممصممومرنوط گم توتلاماحاعصا ,۷۷۰ رتتتطامحط2 0ص .۲ رامعم .. 


یکلا ر5و۲۲۵ و۷۵ نت۱ فاصم ردعامعوع ۷ ۲۱۵00 موتقا ۶ه ومنصهطمع۱ ۳4 مط1. رل .1 و۵016 . 
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و.گ .1 7الحتصه۲۱ مه و.؟ .ظ ات۲ وربا .ظ ۷۷۵۵05 . 


1427, )1969(: 


فا نی چا 


10011060 


ص 


219-233, )1984( 


9 


)1980(, 


سال بیست و پنجم. شماره یک ۱۳۹۳ 


۰ 2 ,۷۲۵61 ۲۱0 .ل و"عصهام 0عصنامصاً صچ موم ما تقصتجصعا ما ممتلهصص۵] ۲۷۷2۵۲۷۵ ربظ .1 مطتصه‌زصه‌ظ 


554-573, )1957(, 


۷ 


7 


.(1963) ,321-334 .00 و(6)3 وکزیا .عاظ ر"عصعام 0مصنامصا مج م0 ما 0تنوتا ۵۶ وتاطاه‌اک رو بطا۷ .8 


۷۵۲ ۱۵۷ ر05010068اض1 رقصمتاهلآنمرد از معط لهمزاممم‌ط/۱2 م1 :۷۷/۵۷۵۵ تماه ۰۷۷ ول .۲ رت801 .9 


(1957) 
۳۵56 الوصا معتطصصعت رد 2صمحمصماظ اهمتووطم۵ع0 0۲60صنک 4ج ۲۷۷2۵۷۵۵ صهعم۵ و۷ مطفتمن 
(1934) ,10800 


ماتتاعصا عاععصیجی ر"قامصصعطل معمای ما ممماهتنامتگومم و۷۵۷ ماصاماو ۵۶ ممتاجعه0۳00 م1" رنه .۲ مققتطمط 1[ . 


1937(۰) ما۲۱ رزعمامصطمع 1 0۶ 


٩۵ ۷۵11, )1938(‏ بألل۷۵0۵۲-۳1 را فتممصلمم متلنتجی ۲۱۲ بو ومتصهطمع]۱۷ 8نبا ۰۳ رویط عنام . 


10۷76 امصصقطه 0مصنامصا صا و۷27۵ لام ۵۴ حصمامامی مطا گم ممتانامی هممتممهمطاه۱۷ ر.۳ بر رتماوومرن1 . 


.(1949) ,149-194 .00 ,2 ,۷۲۵625 .ام ع با( ,ناوجون 


,1 . ,"۲2101211 ۱2اه تملصه ما اعفطه گم ومتممطمع۰۱ ر.ت,۲۱ رب آمممه ۷۷ 220 .۷.۲ رظ۷۵0۵ . 


(1971) ,1367-1386 .00 ,(9769 ,.«:ظ 


۰ ۷011 ۱۱7 ,۷۷۱۱۵ و(تمات0ع) ٩.‏ ۷۰ رانک رع۲۱۷۵۲01۵ ممماوا۲۱۱1 ۶۱0 مصها۵۷۵ .۷۷ ۷۷۰ واامصطط۲ . 


145-175, )1978(, 


ر0۳) کعاا۵ 0۲( ۱۲۵ رعاننلجمه متاقصصکان مممای ظا ۱۷۵۷۵۵ ۵ 0۳۵۵2۵۵0 1" ربق .۲ مفقحطمط[ . 


,(1939) ,214-7229 .وم ,۵۶102 .بنعلنا 


(1970) ,2565-25860 .00 و(96)12 .۲۷ ,۲۷۵۲۵ ,۷۵۷۵۹۲ ۲۵۱۱ کصمصمططیهج ۳۵۳1۵ ریک یک رکام10ظ . 
و(10)2 .۳9 ۷0۳۵0۰ .ل ,10۷7۲ 0صهاو7ن جر م۷۵7۵ (امط ر.ظ ۷۰ رطعصاه 200 وت بط ولا وبا مققطن ور .۵0 وتان . 


00۰ 110-117, )2005(, 


۳۲۵۵۱۵ ]۵ و۷00۱ ۱۷۱۷۸۱۰۱۵۲۱۵۵۱ ۱۳۵6۰ "ومصذامض1 چم 10۵005 صر ۱۷۵۷۵۵ ۰۵11 و.ظ رتعنا0۲0ط0ظ . 


(2007) ,361-366 .00 و(60100156) .ظ مطه/12ظ و.ظ رمتله621012-27) رکععا۳۵5۵0 ۲۸۷۵۲۵۲ 
چاه صرح وامصصفطه صمورن 0مصتامصا مد قم72۲ ۲۵ ۵ تاو اهمزمصصباه م1 ویک وتاک2 200 ۷۰-۲۰ بقتا6 


.(2006) ,1003-1027 .00 و(50)9 وک8تیا ۱۷/6۵۰۰ ۱۷۵۱۸۵۲۰ ۷ ۳۱ و حتاحصنط ۱۷۷۵۷۵ 


یکلا وتعافعطمنط ۷۷۱۱۵ ردوم۲۱ #«مالمطه مممهعع۳ مهد ومطام ممتسطمرمن ماو وی رمتم1 . 
2001(۰) 
جع وهمامصنا تم فممتاتامی لهمزمصصتالط ر.ل .1 روما0ع۳ مضه ر.رت) .۳۲ ۸ 9 رقلله۲ ویک .ظ ۳00 


و396 .۷۲6۵ با .ل عقاو 0620و ج ۵۶ بانازماحاعص ۲۵۱۱-۷۵۷۵ قصح ممتتتاآمن :وم‌هابط ماماتومرعلاهع ص1 ۲۱۵۴۲5 
(1999) ,223-256 00۰ 
۲ 5۳0211077 [1082وصمحصیل-۲۵ مطه متام ۲۵۲ 0مطاعصصر ۸ ۷۷ مطه صرح قصصمامامنم صصحصصمن؟ ر.۲ .۲ روت0 1 


(1992) ,43-68 .00 و338 ی 56۲۰ 10۳۵00۴ 506۰ ۰ 1۳۵۳6۰ .ک90اظ و" عصمتاعناوه 


10 


20 


بم 


2 


22. 


23. 


1-0 ۷۵۷۵ وتحالامد مص‌لمم:ها مصه ممن‌اممره-صمه ۵ عصتام۰۱۵0 ول رتاطا0بوطاماماه 1 صرح بخ رل۱۷]2۴02 .24 


2011(۰) ,945-955 .00 ب(15)6 .و از .لگ بلنص مدز ۰۲۱۷۱۵۰۲۷۷۸۵۳ عصتننامردعا0مطاو عصتونا 


سال بیست و پنجم. شماره یک» ۱۳۹۲ نشریه مهندسی عمران فردوسی 


۸ کاربرد الکوی ۷۷۸ در شبیه‌سازی عددی امواج غلتشی.. 


وله 4ج فصممناوم مماه7 ۲مالمی عصتویا ممتاما ای ومع تام ریخ رلاته‌طماصما همه .ظ رطع001ع2 .25 

2002(۰) ,610-622 00۰ و(40)5 ,.5ع1 1۳۷۵۲۵۷۱۱6 ۰ ,"00طامحظ ناگ م2۷۵۲ 
فنص صرح مصصنام۲ عاتص موق گم توقای ممواهم‌تتمن ولا رقلطموع]. 24 ,۷۲ ر1/016۵60۷78-1۷]607100۷7۵ .26 
۲ و0 1010 مه ومجنهووه) مها حز ۱ کممتهناوع ماج هجو فص ۶۵ 05مطامجط هام 


:(2006 ,874-891 .00 ,(86)11 ,-/1۷۵۵6 ۱0۵۰ 0عزاوهچق . 


نشریه مهندسی عمران فردوسی سال بیست و پنجم. شماره یک ۱۳۹۲ 


